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We study the properties of QCD at high baryon density in a finite volume where color super- 
conductivity occurs. We derive exact sum rules for complex eigenvalues of the Dirac operator at 
finite chemical potential, and show that the Dirac spectrum is directly related to the color super- 
conducting gap A. Also, we find a characteristic signature of color superconductivity: an X-shaped 
spectrum of partition function zeros in the complex quark mass plane near the origin, reflecting 
the Z{2)l X Z{2)r symmetry of the diquark pairing. Our results are universal in the domain 
^ L <C where L is the linear size of the system and is the pion mass at high density. 



PACS numbers: 12.38.Aw, 21.65.Qr 

Revealing Quantum Chromodynamics (QCD) in the 
regime of finite temperature (T) and chemical potential 
(/z) is important for understanding a wide range of phe- 
nomena from ultrarelativistic heavy ion collisions, the 
early Universe, and neutron stars to possible quark stars 
[H. A lot of theoretical progress has been made by the 
first-principles lattice QCD Monte Carlo simulations in 
the study of the finite-T regime Q. However, the appli- 
cation of the lattice technique to QCD at finite /i is still 
hampered by the notorious fermion sign problem: cal- 
culation of the QCD partition function requires dealing 
with a path integral with a measure including a com- 
plex fermion determinant. This is the main reason why 
our understanding of the properties of QCD at finite fi 
is still immature, except at asymptotic high /i where the 
ground state is shown to be the most symmetric three- 
flavor {Nf = 3) color superconductivity (CSC), i.e., the 
color-flavor locked (CFL) phase 0, 0] by using the weak 
QCD couphng calculations. 

In this Letter, we demonstrate exact analytical results 
for QCD at high /z specific for a large but finite volume. 
By matching the partition function of QCD against that 
of the effective theory of CSC, we derive exact sum rules 
for the Dirac eigenvalues (Dirac spectrum) as well as the 
spectrum of partition function zeros (Lee- Yang zeros [H|) 
in the complex quark mass m-plane. As is well known at 
/i = 0, the Dirac spectrum and the Lee- Yang zeros spec- 
trum of QCD in a finite volume are very closely related 
to the chiral symmetry breaking: the chiral condensate 
(qq) is directly connected to the Dirac spectrum by the 
exact relations such as the Banks-Casher relation [6| and 
the Leutwyler-Smilga sum rules [7.]. Also, a nonzero (qq) 
implies the existence of a line of the Lee- Yang zeros going 
through m = in the complex m-plane independent of 
H Q. Nevertheless, such exact relations at finite /i and 
the relevance of the Lee- Yang zeros spectrum to the CSC 
have not been fully understood. 

As we shall show below, the Dirac spectrum at high 
fi is intimately related to the CSC gap A, rather than 
to (qq), through our spectral sum rules. Also, a nonzero 



gap A necessitates an X-shaped cut of Lee- Yang zeros 
near m — 0. In particular, the Z(2)l x Z{2)ii symmetry 
of the diquark pairing plays a crucial role on both Dirac 
and Lee- Yang zeros spectra. Together with the exact 
results at /X = [H, we expect that our results impose 
strong constraints on their possible spectra and provide 
important insights to the properties of QCD at finite fi. 

In the following, we will focus on QCD with Nf — 
3 (light up, down and strange quarks) at finite quark 
chemical potential fj. living on the four-dimensional torus 
Vi — LxLxLx/3 with /? = 1/T ^ L. Let us consider 
the Euclidean QCD Lagrangian defined as £qcd = 9(2?+ 
M)q + Cg with Cg = jF^^F^" and the Dirac operator 

V = j^id^ + igA^) + /i7o. Here q is the quark field 
and Af^ = A^^t"" is the gluon field with the color SU{3)c 
generators t° (a = 1, 2, • • • , 8). M is the complex quark 
mass matrix, g is the QCD coupling constant and F^^ — 
dfiAi, — duA^ + ig[Af^, A^]. Since T> is not anti-Hermite 
with /i > 0, its eigenvalues iA„ are generally complex 
values, whereas iA„ are pure imaginary at /i = 0. Even 
so, 2? preserves the chirality, {75,2?} = 0. The chirality 
ensures that if ?A„ is the eigenvalue of V {'D^n = i^nipn), 
then so is — iA„. 

The QCD partition function Zqcd involves a sum over 
the different topological sectors of the gauge-field config- 
urations characterized by the integer topological charge 1/ 
as Zqcd = Y^^e'^^Z^. At ^ > Aqcd (Aqcd: the typi- 
cal scale of QCD), however, the topological susceptibility 
is highly suppressed as {v'^) cx (Aqcd/m)^ owing to 
the screening of instantons in the medium together with 
the asymptotic freedom of QCD. Thus, we can focus on 
the topological sector 1/ = alone. 

The QCD partition function with h' — can be writ- 
ten in the functional integral using the symmetry iA„ <-!■ 
—iXn- 
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is the average of O over all gauge configurations with 
Sg being the classical action of the gluon field. Zqqd is 
normalized so that Zqqd — 1 when quark masses are 
turned off. 

We shall give the partition function Z-eft from the 
effective theory of the color superconductivity (CSC). 
For definiteness, we consider the most predominant di- 
quark pairing, the color- flavor locked (CFL) phase [3]: 
{{<lL)lC{qLf^) ~ eabct^ok[dlU and {{qB)iC{qR)l) - 



where A: {a,b,c) are the flavor (color) 



indices, and C is the charge conjugation matrix. The re- 
markable feature here is that chiral symmetry is dynam- 
ically broken by the diquark condensate: the symmetry 
breaking pattern of the CFL phase at asymptotic high 
^i is SU{3)c X SU{3)l x SU{3)r x U{1)b x U{1)a ^ 
SU{3)c+L+B. X Z{2)l X Z{2)r [The Z{2)l x Z{2)r sym- 
metry left reflects the fact that we can change the sign 
of the left-handed or right-handed quark fields indepen- 
dently]. As a result, we have 8-1-1-1-1 Nambu-Goldstone 
(NG) modes associated with the breaking of chiral sym- 
metry, U{1)a and U{1)b symmetries, which we will re- 
fer to as pions, 77', and -ff, respectively. In the follow- 
ing, we will not consider H since its dynamics decou- 
ples. In the CFL phase, we also have gluons and quarks; 
the gluons acquire a mass comparable to the CSC gap 
A (X [m, [13] by the Anderson-Higgs mechanism 

when the SU{3)c symmetry is broken; the octet (singlet) 
quarks of the unbroken SU{3)c+l+b. symmetry have the 
mass gap A (2A) 3. 

We then specify the microscopic domain (or e-domain) 
of the CFL phase. The microscopic domain of QCD at 
/i = is specified by A~^ L <^ "^n^i where mn is the 
mass of pions at fi = and A is the mass scale of the 
lightest non-NG modes (i.e., the p meson mass nip) 
The corresponding microscopic domain of the CFL can 
be defined as 



The CFL effective Lagrangian in the Minkowski space- 
time up to the leading order 0{M^) is given by [14| 



Ceft = ^Tr[VoSVoSt - i;2a^sast] 
■ 'f- [doVdoV* - vl,d,Vd,V*] 
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T/(TrMS^)2 - T/Tr(MEt)2 ^ ^ 



where S = exp(i7r''A''/A) and V = eyi])[2irj' / {^f^,)] 
are the pion and rj' fields respectively, /tt [f-q') is the 
pion {rj') decay constant, (^^r;') is the pion {q') ve- 
locity, A° (a = 1, 2, • • • , 8) are the Gell-Mann matrices, 
and the covariant derivative including the effective chem- 
ical potential (Bedaque-Schafer term 15|) is given by 

VoS = dot + i (^) t-it (^) with the Fermi 

momentum pp- The quantities /,r and can be pertur- 
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batively computed at ^ Aqcd as 



Pf 



In Eq. ([3]), we have neglected the mass term of or- 
der 0{M), since this term originates from the instan- 
ton contribution and is suppressed at asymptotic high 
/i [10]. Thus, the leading mass term in the CFL effec- 
tive Lagrangian is ©(M^), unlike the 0{M) term in the 
usual chiral Lagrangian at low /i. A more intuitive ex- 
planation for this fact is that Mqq is prohibited by the 
Z{2)l X Z{2)r symmetry, but {MqqY is not. 

In the domain L <^ one can neglect the contri- 

bution of the kinetic term. Then the partition function 
for the CFL effective Lagrangian reads: 

f { 3A2 

Zeft= /rfSexpfy4— [(TrMEt)2_Tr(Afl]t)2] detS-|-H.. 
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where ttItt is the mass of pions associated with the CSC at 
high [I. The first condition in Eq. ^ follows by compar- 
ing the contribution to Zeft of the pions, e~™'"^, to that 
of the other heavier particles, e^^^ . This condition al- 
lows us only to deal with the pions described by the CFL 
effective Lagrangian. On the other hand, the second con- 
dition in Eq. ([2]) means that the Compton wavelength of 
the pions is much larger than the linear size of the box, 
so that the CFL effective Lagrangian can be truncated to 
its zero momentum sector. Note that the second condi- 
tion is automatically satisfied at sufficiently high [i with 
h and quark mass m fixed, since m-,^ Am/ 11 (see Eq. 
^ below) together with the relation A '-^ /iexp (^—^^^ 

[l3]. Note also that, in the domain ([2]), temperature T is 
low enough for the CSC to be realized since T A ^ Tc 
with Tc being the critical temperature of the CSC. 



where the integral is over E = 'EV £ U{3) and Zeft 
is normalized so that Zeft = Zqcb = 1 in the chiral 
limit. In Eq. we have also neglected the effect of the 
effective chemical potential. [If one includes it, Eq. ^ 
can be expanded in terms of not only (yA'^)'^0{M'^) but 
also VO{M'^). In the domain A^^ ^ L, however, the 
latter is negligible.] 

Owing to the property of E e C/(3), [(TrMEt)2 „ 
Tr(MEt)2]detE = 2detAf Tr(M-iE), Eq. ® can be ex- 
pressed analytically as shown in Ref. [16|. In particular, 
in the flavor symmetric case hi — ml, one can evaluate 
Zeft in a simpler form using Weyl's formula [tI. Ii7|: 



Zeft = det [Ij-i{x)\ 

Q<i,j<2 



(5) 



where Ii,{x) is the modified Bessel function and x = 
ZViVi? A^ / TT^ . The combination of m2A2 is expected, 
since m?' acts as a source for A2. It should be remarked 
that the expression © is exactly the same form as the 
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partition function with = at /i = which is 
given by Eq. ([5|) with the replacement of the argument: 
X ^ x' ~ V4m|(5q)|. This is a novel correspondence be- 
tween the CSC phase and the hadronic phase, and may 
have relevance to the idea of their continuity [l^. In 
Table. [H we summarize our main results at /i ^ Aqcd 
below compared with the results at = in Ref. 01 • 

Expanding in terms of quark mass M and performing 
the group integral over S G J7(3) order by order, Eq. (jU 
reduces to the following form up to ©(M^): 

^EFT-1 + \ [y^^ [(TrMtAf)2-Tr(MtAfAftM)].(6) 



Using the relation, det[l+e] = l+Tre+i[(Tre)2 - Tre2] + 
©(e"^), one can expand the QCD partition function ^ in 
terms of the quark mass matrix M . Then one obtains 
the spectral sum rules for the Dirac eigenvalues zA„ by 
matching this expansion against Eq. ([6]). By rescaling 
Zn = -y/V^AAn, the results read 





(7) 



(8) 



where the summation is taken over z„ satisfying 
Re(z„) > and |z„| < y^A^ (|A„| < A). These rela- 
tions are highly nontrivial, since sums of inverse powers 
of complex A„ with the average taken over the gauge con- 
figurations give the real value involved with the CSC gap 
A. In particular, the sums in Eq. ([8]) are identically zero, 
which is a direct consequence of the Z{2)l x Z{2)r sym- 
metry of the quarks at high /i. This situation should be 
compared with QCD at /i = (at low fi [19]), where 
sums of inverse powers of real A„ (complex A„) take pos- 
itive values involved with the chiral condensate (qq). As 
in the case oi fi — 2^, 21 1 , our spectral sum rules must 
be universal (i.e., independent of microscopic details) in 
the domain ([2]). 

By using the spectral density p{X) — ((X^n ^'^{^~' ^n))), 
the first sum rule in Eq. ([7| reduces to the relation: 

where C+ = {z : Re(2:) > 

This implies the existence 
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0} and (Pz — d(Kez)d(lmz 
of the microscopic limit of the spectral density defined as 



Psiz) 



lim 



V4A 



(9) 



From Eq. ([9]), we find that the microscopic spectral den- 
sity at high /i is governed by the CSC gap A. Also 
Eq. ^ shows that the linear spacing of eigenvalues (5A„ 





order parameter 


source 


SXn 


Lee- Yang zeros 


M > Aqcd 






oc I/VV4 


Re(m) — ±Im(m) 


/i = 


(qLqR) 


m 


oc I/V4, 


Re(m) = 



TABLE I: The summary of our main results at /i S> Aqcd 
compared with the results at /i = in Ref. 0]: order param- 
eters, the corresponding source terms, the linear spacing of 
eigenvalues SXn, and the spectra of partition function zeros 
(Lee- Yang zeros) in the complex m-plane, respectively. 



in the complex A-plane satisfies SXn oc l/^/V^. Since 
dXn oc I/V4 at low fi with (qq) ^0 0, and (5A„ oc l/V^^"^ 
in a free theory, our result indicates a sizable deformation 
of the Dirac spectrum due to the dynamics of the CSC. 
We expect that the random matrix theory (RMT) [2l| 
or the supersymmetric approach [2^ incorporating the 
CSC and the symmetries of the CFL not only reproduce 
the above results, but also clarify the concrete form of 
the microscopic spectral density ps- 

Let us consider the partition function zeros (Lee- Yang 
zeros) in the complex m-plane in the flavor-symmetric 
case. Using the asymptotic form of I^{x), we find the 
partition function ([5]) for |a;| ^ 1 as 



^EFT 



{x 



-5/2 



cos 4z 



Stt 

T 



(10) 



Then, the Lee- Yang zeros for |a;| ^ 1 are given by a: = 
-i(n -I- |)f (n g Z). Remembering x = 314?Ti^A^/7r^, 
the zeros are spaced along the lines 



Rc(to) = ±Im(TO), 



(11) 



at /i 3> Aqcd- In the thermodynamic limit V4 00, the 
density of the zeros increases and they join into a cut in 
the vicinity of massless limit m = 0. However, the cut 
does not go through m = 0, since dx/dm = 6V4mA^/7r^ 
and the density of the zeros vanishes at m = 0. 

In Fig. [U we draw spectra of the Lee- Yang zeros in 
the complex m-plane near m = 0: the spectrum (a) at 
fi ^ Aqcd is the exact result obtained above. For com- 
parison, we show the spectra (b) ai p, — previously 
obtained exactly [2^, and (c) for fi > iic obtained from 
the RMT 0, where iic is the critical chemical potential 
of chiral symmetry restoration. In case (b), the density 
of the zeros at m = is finite as dx' /dm = Vi\{qq)\, 
and the chiral condensate {qq) takes nonzero value. As 
/i increases, the result (c) shows that the zeros will move 
away from the origin and chiral symmetry restores 

Our exact result, however, demonstrates that the sce- 
nario (c) of the RMT should suffer from dramatic mod- 
ifications if the effects of the CSC is taken into ac- 
count: even at /i ^ Aqcd, there exists an X-shaped 
cut (fTTj) near m = as shown in Fig. [TJa), reflecting the 
Z(2)l X Z(2)_R symmetry of the diquark pairing. The 
distinctive feature of the cut (fTTj) is that there is a dis- 
continuity in the m'^-plane at m = along the real axis 
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but not in the m-plane, which is related to the fact that 
A ^ and {qq) = at /Lt Aqcd ■ How the spectrum of 
the Lee- Yang zeros evolve from (b) to (a) as increases 
is beyond the scope of this paper and needs further inves- 
tigations. Note that summing over the topological sector 
J/ would also change the scenario (c) where i/ = alone 
is considered. 





Re(ffj) 



(a) 



lm{m) 



tlm(?w) 



Ke(in) 



Re(m) 



(b) (c) 

FIG. 1: The spectra of the Lee- Yang zeros in the complex 
m-plane in the vicinity of m = in the thermodynamic limit 
14 — > CX3: (a) A* ^ Aqcd (exactly obtained in our analysis), 
(b) ^ = ^3], and (c) fi > fj^c with the critical fic of chiral 
symmetry restoration obtained from the random matrix the- 
ory The density of the zeros fades away as m in case 
(a), while that remains constant in case (b). 

It is important to generalize our spectral sum rules 
or to directly investigate the distributions of the parti- 
tion function zeros at lower baryon densities. One can, 
e.g., match the QCD partition function at finite density 
against the effective theory of the generalized pions [24] 
in the entire span of the density where the microscopic 
regime can be defined as nip^ <C L ^ according 
to Ref. ^li*!. Also the generalization of our spectral sum 
rules to QCD-like theories, such as the two-color QCD at 
high density, would be an interesting problem to be in- 
vestigated 1255, which can be tested on the lattice QCD 
simulation. 
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